Astronomy & Astrophysics manuscript no. 2279manu 


© ESO 2009 


May 27, 2009 





Letter to the Editor 



o 
o 

(N 
>> 



(N 
< 

a 
^' 

Oh; 
O 



> 

in 
m 

o 

On 
O 



O and Na abundance patterns in open clusters of the Galactic disk 

G. M. De Silva', B. K. Gibson^, J. Lattanzio^ and M. Asplund"^ 
' European Southern Observatory, Karl-Schwarzschild-Str 2, D-85748 Garching, Germany 

- Jeremiah Horrocks Institute for Astrophysics & Supercomputing, University of Central Lancashire, Preston, PRl 2HE, UK 
^ Centre for Stellar and Planetary Astrophysics, Monash University, Victoria 3800, Australia 
* Max Planck Institute fiir Astrophysik, Postfach 1317, D-85741 Garching, Germany 



Received ; accepted 



ABSTRACT 



Aims. A global O-Na abundance anti-correlation is observed in globular clusters, which is not present in the Galactic field population. 
Open clusters are thought to be chemically homogeneous internally. We aim to explore the O and Na abundance pattern among the 
open cluster population of the Galactic disk. 

Methods. We combine open cluster abundance ratios of O and Na from high-resolution spectroscopic studies in the literature and 
normalize them to a common solar scale. We compare the open cluster abundances against the globular clusters and disk field. 
Results. We find that the different environments show different abundance patterns. The open clusters do not show the O-Na anti- 
correlation at the extreme O-depletion / Na-enhancement as observed in globular clusters. Furthermore, the high Na abundances in 
open clusters do not match the disk field stars. If real, it may be suggesting that the dissolution of present-day open clusters is not a 
significant contribution to building the Galactic disk. Large-scale homogeneous studies of clusters and field will further confirm the 
reality of the Na enhancement. 
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1. Introduction 

The O-Na abundance anti-correlation observed in Galactic glob- 
ular clusters (GCs) is tho ught to be an intrin sic property of 
the clusters (see review bv lGratton etani2004 and references 
therein). It is considered a global feature of GCs, because the 
anti-correlation has been observed in all GCs subject to high- 
resolution studies (Carretta et al. 2006) and is present among the 
cluster giants, as well as unevolved stars dCarretta et al] 120041: 



iRamfrez & Cohenll2002t lOratton et al.ll200C Iron and the heav 
ier eleme nts show little star- to- star abundance variations within 
a cluster ( ISuntzeftlll993 l). The currently supported theory for the 
origin of the O-Na anti-c orrelation is that of primor dial pollu- 
tion from previous stars dCottrell & Da Costal 1 1 98 1[) . although 
the mechanism responsible for the pollution is still unclear. 
Several hypotheses are presently bein g debated, including pol- 
lution by the ejecta of massive stars jPecressin et al. 2007a'' 



or by AGB stars undergoing hot bottom burn ing dVentura et al 
IMulFeimerern]l200llKarakas et alJl2006h . 

In general, the halo stars have a chemical composition sim- 
ilar to that of the GC s with the exceptio n of the lighter ele- 
ment abundance trends jGratton et alj2000 ). That the O-Na anti- 
correlation is not seen in the halo stars presumably reflects the 
different chemical evolution of the high-density cluster environ- 
ment. Therefore stellar populations showing GC-like chemical 
evolution can be isolated from ot her popula tions by examin- 
ing the O-Na abundance pattern CGeisler et al]|2007l) . The abun- 
dance mismatch with the halo has important implications for the 
formation of the Galactic halo. It is possible the halo was built 
up either from clusters that have long since dissolved, which had 
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a different chemical enrichment history to the surviving GCs, or 
from the stars lost from the present-day GCs before taking on 
the effects of the polluting stars. 

A similar investigation can be applied to the progenitors of 
the Galactic disk. It is often stated that open clusters (OCs) are 
the object s of choice for tracing the star-formation history in the 
disk (e.g. lFrielll995h . Present-day OCs cover a wide range in 
age, metallicity, and position. They are used to examine the disk 
chemical evolution, from the solar neighborhood to the outer 
Galactic disk. Older OCs are considered particularly useful be- 
cause they provide a time line for change. In general the OC 
abundances match that of the field, a lbeit with a larg er scatter 
(iFriel et al.ll2002t iDe Silva et al.ll2007l: ISragaglia et alJl2008l). 

It is interesting that some OCs are Na-enhanc ed dFriel et al] 

120051: iJacobson etaT] |2007: Brag aglia et ai]|2008l) . Whether this 
is truly an intrinsic property of the clusters or an artifact of 
the abundance measurements remains unclear in the literature. 
A discrepancy in abundance patterns between the OCs and the 
disk has important implications for our understanding of disk 
formation. Recall that the O and Na abundance anomaly is the 
characteristic difference between the GCs and the halo. To our 
knowledge no such characteristic abundance patterns have been 
reported for OCs. In this Letter we explore the abundance pat- 
terns of O and Na among the OC population and compare them 
against the GC abundance anti-correlation and the disk field 
abundances. 



2. Data samples 

We combine high-resolution studies that derive both O and 
Na abundances in open clusters. The homogenized open clus- 
ter abundances, associated rms scatter, the number of stars 
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per cluster, and references are given in Table [T] To represent 
the GC O-Na anti-correlation, we use abundances of red gi- 
ants in NG C 2808 as a templat e of the global GC O-Na anti- 
correlation dCarretta et alJ 120061 see Fig. 5 ). For comparisons 
with th e disk field, we use the compilation bv^S oubiran & GirardI 
(l2005h mainly base d on fie ld dwarfs, the abundances pubUshed 
by iMishenina et af] ((2006) based on field red clump stars, and 
the disk giants studied by .Fulbright et al.. (.2007,) . 

Table 1. Open cluster sample 



ID [Fe/H] N [0/Fe] cr [Na/Fe] cr Ref. 

T" 
2 
2 
2 
3 
4 
5 
6 
6 
2 
5 
4 
7 
7 
8 
8 
10 
4 



References: (1) JFriel et alj t200^; (2) lYong et alj HMB): (3) ICarretta et alj t200^; (4) 
IPace et aljr2008t): (5)IGratton & Contarinilll994l): (eterown et aljti996l): (7)ICarTetta et al" 
<2007hl) : (8) IJacobson et alj 120071) : (9) iJacobsonetalJ I2O08C - (10) ITautvaisieneTTar 

As always when comparing different literature sources, the 
presence of systematic effects must be highlighted. As the sys- 
tematics arise at various stages, it is very cumbersome to ac- 
curately correct for all possible effects. For example, the stellar 
effective temper atures (Teff) could be subj ect to large system- 
atic difi'erences. ICarretta et all (l2007bl) and lBrown et al.1 (1 19961) 
derive T^ff based on photometry, while other studies are based 
on spectroscopy. A difference of lOOK in Tetj may result in Na 
abundance differences up to 0.1 dex. Since the studies are on 
different clusters, such systematics are difficult to quantify. To 
homogenize the different studies, we have normalized the pub- 
lished abundances to the new solar values of loge (O) = 8.70, 
loge (Na) = 6.21 and loge (Fe) = 7.51 (Asplund et al. 2009). 
This normalization was applied to the abundances of all OCs, 
GCs, and field stars, except for studies based on a strictly differ- 
ential analysis with respect to the Sun. 

Most of the studies use the same elemental lines to derive 
oxygen and sodium abundances. For oxygen, either one or both 
of the [O I] forbidden lines at 6300.3 A an d 6363.8 A ha ve been 
used, where NLTE effects are neglig ible (lAsplundll2005l) . While 
all othe r studies employ ful l spectral synthesis to deri ve O abun- 
dances, ICarretta et al] ("2006) and iPace et al.1 (|2008|) use EWs, 
taking the effects of line blending into account. 

For sodium, most studies use Na I lines at 5682.6, 5688.2, 
6154.2 and 6160.7 A. The gf values used vary among the stud- 
ies. In many cases this difference is small with negligible effects 
on a bundances. The greatest differences are by Jacobs on et al.l 
(12007. Be 17 & NGC 6939), where the adopted log gf values 
are higher by 0.2 - 0.3. As this results in a larger Na abundance 
of about 0.2 - 0.3 dex, we adjust their published abundances to 
account for the offset in the gf values. 

For the open cluster data, we use only published Na abun- 
dances based on LTE analysis in order to minimize systemat- 
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Fig. 1. Normalized [Na/Fe] vs. [O/Fe] for open clusters (red cir- 
cles) compared to globular cluster NGC 2808 (green squares). 
The red triangles are open clusters based on solar-type stars by 
iPace et all (l2008l). The black line rep resent a Galactic chemical 
evolution model dHughes et alj|2008t) . 



ics. The Na lines are, however, subject to NLTE effects, where 
the NLTE corrections are dependent on the stellar parameters , 
[Fe /H], and the Na abundance (e.g. iMashonkina et al.l 120001 : 
Tak eda et alj |2003). While the cluster sample spans almost one 
dex in [Fe/H], we do not see an Na abundance trend with metal- 
licity when assuming LTE. Such a trend would be expected in 
more metal-poor stars where the NLTE corrections are larger. 
For the globular cluster NGC 2808, the authors correct their 
Na abundances for NLTE effects following Grat ton et al. (199^. 
The Na abundances for the field clump stars by IMishenina et all 
(2006) is based on NLTE analysis. In all other cases we adopt 
LTE abundances. 



3. Results and discussion 

In Figure [1] we plot the normalized [Na/Fe] and [O/Fe] for the 
sample of OCs and GCs. Overplotted is the Galactic chemical 
evolution (GCE) m odel that corresponds to the one described by 
Hu ghes et al.1 ("2008). Only the final ~4Gyrs (for [Fe/H] > -1) is 
plotted for the annulus corresponding to the solar neighborhood. 
The plotted values for OCs are the mean cluster abundances, 
whereas for the GCs the individual stars are plotted. We chose 
to plot only the mean value for each open cluster because no 
significant sta r-to-star abundances variations are found within an 
open cluster ( De Silva et al.ll2006l) . 

The OCs sit in the middle range of the O-Na anti-correlation 
of the GCs. They do not reach the extreme values of O-depletion 
/ Na-enhancement and vice versa. If we disregard the only two 
clusters with subsolar Na abundance (Praesepe and M67), this 
might look as if the OCs follow the GC anti-correlation, where 
all stars in a given OC have abundances similar to an individ- 
ual star in a GC. The origins of the anti-correlation in GCs is 
thought to stem from primordial pollution. A plausible source 
is intermediate mass AGE stars undergoing hot bottom burn- 
ing where the temperatures are hot enough to activate the CNO 
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Fig. 2. Surface gravity v s. Teff for open cluste r stars com- 
pared to disk giants from iMishenina et alj (l2006h (green) and 
iFulbright et alj (l2007h (blue). Filled red circles are open cluster 
stars with T^ff > 4500 K and log g > 2 and the open red circles 
are the remaining open cluster stars. 



and NeNa cycles dKarakas & Lattanzioll2003l) . Here Ne is used 
to produce Na via proton-capt ure, while O is depleted in the 
ON cycle dKarakas et al.1 12006'). These AGB stars then pollute 
the proto-GC via low-speed winds, which are retained in the 
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Fig. 3. [Na/Fe] and [O/Fe/] vs. [Fe/H] for open clusters and 
field. Red filled circles are cluster abundances from stars 
with Teff > 4500 K and log g > 2. Red tr i angles are clus- 
ters based on solar type stars by iPace et all d2008l) . and the 
red open circl es are a ll other open clusters. Overplotted are 
disk dwa rfs bvlSoubira n & Girard (2005) (crosses), disk clump 
stars by Mishenina et alj (I2OO60 (green), and disk giants by 
iFuibright et al.. C2007.) (blue). 



cluster to form new stars that are O-de pleted and Na-enhanced 
dVentura et al.ll2002l:lFenner et al.ll2004D . 

Since all stars are chemically homogeneous in an OC, 
they must have form ed from the same (well mixed) gas cloud 
dPe Silva et al.l2066l) : therefore, if pollution from an earlier gen- 
eration has taken place, it has to pollute all stars homogeneously 
or take place before the cluster assembled. It is possible that a 
few high-mass stars form shortly after the cloud assem bles and 
enrich the cloud fairly uniformly dMcKee & Tanll2002h . In such 
a scenario the pollution effects are observed in the proceeding 
generation of stars homogeneously within the cluster 

Another possibility is that the first generation of OC stars 
were dispersed and only the polluted stars are presently bound 
as a cluster. The level of pollution and the degree of O-depletion 
and Na-enhancement of the polluting stars must vary at the dif- 
ferent site and time of the formation of the OCs in order to ac- 
count for the abund ance spread acro s s the O C sa mple. 

The studies of iBragagha et"!!] d2008l) and iDe Silva et all 

d2007h also suggest that Ba could be enhanced in OCs. The real- 
ity of this enhancement requires further investigation, but it may 
indicate other sources resp onsible for pollutio n in OCs. With the 
exception of NGC 1851 dYong & Grundah]||2008l) . abundance 
variations of s-process elements in GCs are found to be in- 
significant and unrelated to the O-Na anti-correlation. The likely 
source of s-process elements are low mass AGB stars, whereas 
more massive AGB stars are thought to be responsible f or the O- 
Na an ti-correlation in GCs through hot bottom burning dRenzinil 
120081). 

Figure [T] could be interpreted as the OC population showing 
an O-Na abundance correlation, albeit with a wide spread. The 
GCs also show a broad spread about the mean anti-coiTelation. 
Fitting a least square regression line provides a coiTelation co- 
efficient of 0.5 when all OCs are considered, or ~0.78 when 
the data is split into two groups with a vertical offset of 0.3 
dex. The GCE model also shows an abund ance correlation. Th e 
model corresponds to the one described bv Hughes et alj d2008l) . 
A 'dual infair (halo + disk) framework for the solar neighbor- 
hood was employed in which the infall of primordial gas during 
the halo phase occurred on a rapid ~50 Myr timescale; a second 
(disk) phase was delayed by ~1 Gyr with respect to the first, and 
then gas enriched to 10% solar meallicity was assumed to infall 
on an ~ 10 Gyr exponential timescale, thereafter. A conservative 
Schmidt Law for star formation was assumed with an efficiency 
constr ained by the present-day gas fraction. A iKroupa et al.l 
dT99l) IMF was employed with lower and upper mass limits 
of 0.08 and 60Mo, respectivel y. Most important, t he ad opted 
stella r yields include those of IWooslev & Weaveii ( 1995L >10 



Karakas & Lattanziol (12003', <8 Mq), and Nomot o et"a!] 



dl997l Type la SNe). We would expect a match in the abun- 
dance trend between the OCs and the GCE model under the 
widely held assumption that the OCs were a major contributor to 
building up the Galactic disk. Slight variations in the model pa- 
rameters would not be sufficient t o recover an O-Na abundance 
anti-correlation dFenner et al.ll2004.) . and recovering such articu - 
lations would require extreme conditions dMarcohni et alj2009l) . 

We now compare OCs against the disk field. To examine 
the stellar parameters of the different samples. Fig. |2] plots the 
log g vs. Tss of the OC stars along with the field stars by 
IMishenina et alj d2006l) and lFulbright et alj d2007h . It shows the 
open clus ter stars with Teff > 45 00 K and log g > 2, which 
match the IMishenina et al] d2006l) stellar parameters. The open 
ci rcles show all other OC sta rs. Omitted from Fig.|2]are the stars 
of ISoubiran & Gira rd (2005), which are mostly dwarfs, and the 
OCs studied bv lPace et al.. (.2008') . which are based on solar-type 
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Stars (Jeff > 5500 K and log g > 4). In Fig. [3] we plot [Na/Fe] 
and [O/Fe] vs. [Fe/H] for the open clusters and disk field stars 
(cf. Fig 4. in Bragaglia et al. 200^. 

The OC oxygen abundances match the various field sam- 
ples, with the exception of Praesepe and Mel 71, which are 
among the most O-depleted clusters. Sodium is significantly 
enhanced in OCs compared to the field samples. As men- 
tioned in Sect. [2l we hav e only used OC LTE abundance. The 
ISoubiran & GirardI ( l2005h compilation also uses LTE analyzes, 
but based on dwarf stars, whereas the OC stars are mostly red 
giant and clump stars. It is possible that the OC giants are show- 
ing effects of internal mixing, m aking them Na-enhanced com- 
pared to dwarfs. Pasquini et a l.l f2004) find Na enhanced by 0.2 
dex in giants compar ed to dwarf members of OC IC 4561, while 
iRandich et all (l2006l) do not see abundance variations between 
main sequen ce and sub-giant stars in M 67. Therefore discrep- 
ancy with the lSoubiran & GirardI (l2005l) sample could be due to 
the choice of tar gets. 

The study by lMishenina et al] (l2006t) is of field clump stars, 
which match some of the OC stellar parameters (see Fig. |2|; 
however, they perform an NLTE analysis using the profiles and 
EWs of the Na lines. They suggest the NLTE correction is about 
0.1 to 0.15 where the LTE values are higher than their NLTE 
equivalent. Further [Sestito et akl (l2008l) find OC Na abundances 
based on NLTE computations produce values lower by 0.1 - 
0.2 dex than using LTE assumptions. Therefore the discrepancy 
between the OC and field clump stars may come from NLTE 
effects. On the other hand, the abundances of disk giants by 
iFulbright et al.l (l2007l) are from an LTE analysis. Therefore the 
Na abundance mismatch between the OCs and this sample is un- 
clear, since here we are comparing similar type stars (see Fig.|2]l 
and both are based on LTE analysis. 

It is important to confirm this possible OC Na enhancement, 
because it has major implications for the formation scenario of 
the Galactic disk. Likewise with the GCs and the Galactic halo, it 
may imply that the OCs present today were not the major build- 
ing blocks of the Galactic disk. Rather, now dispersed clusters 
had a different chemical evolution than those that have survived 
to date. The OC abundances for other heavier elements match 
those of the disk fie ld as is the case between the GCs and the 
halo (ISuntzefflll993 l). The differences between the halo and GCs 
are highlighted mainly in the lighter element abundances, no- 
tably the O-Na anti-coiTelation. 

The presented sample of OCs in this study are mostly old, 
over 1 Gyr in age. Detailed abundance analysis of younger OCs 
are lacking in the literature. It is expected that mos t clusters will 
dissolve into the field within the first billion years jLamers et alj 
I2005h . The survival of extremely old OCs around 10 Gyr in age 
is very interesting as important fossils of the conditions of an 
earlier era. One can question whether such surviving OCs are 
representative of the chemical evolution of the disk, or rather 
whether embedded clusters subject to high infant mortality rates 
are the major contributors to the building of the Galactic disk 
(iLada & Lada 2003). The older, still bound clusters, which have 
deeper potentials, may in fact be the anomalous objects that 
faced a different formation and chemical evolution. 

4. Conclusions 

We have compiled Uterature studies of O and Na in OCs and 
compared the abundance patterns with the GCs and disk field 
stars. We homogenized the various studies to put them on the 
same abundance scale as best as possible. The OC population sits 
in the middle of the well-established GC anti-correlation, and it 



is unclear whether they show an O-Na correlation or a GC-like 
anti-correlation without reaching the extremes of O-depletion 
and Na-enhancement. 

When compared to field abundances, the O abundances of 
OCs match the field stars well. Na is enhanced relative to the 
field. Much of the Na abundance discrepancy could be explained 
by possible internal mixing in the giant stars and NLTE ef- 
fects; however, a possible intrinsic abundance mismatch cannot 
be ruled out due to the mismatch between LTE o pen cluster giant 
abundances and the LTE study of field giants by IFulbright et al.l 
(120071) . Homogeneous abundance studies of clusters and field is 
needed for confirmation. If the Na-enhancement is indeed a true 
feature of old open clusters, it has major implications for un- 
derstanding the formation of the Galactic disk. Further study of 
O-poor clusters is encouraged to determine the O-Na abundance 
behavior at the extreme values of O-depletion. 
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